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Abstract. The formula for dark energy density derived by Gurzadyan and Xue is 
the only formula which provides (without a free parameter) a value for dark energy 
density in remarkable agreement with current cosmological datasets, unlike numerous 
phcnomenological scenarios where the corresponding value is postulated. This formula 
suggests the possibility of variation of physical constants such as the speed of light 
and the gravitational constant. Considering several cosmological models based on 
that formula and deriving the cosmological equations for each case, we show that, 
in all models source terms appear in the continuity equation. So, on the one hand, 
GX models make up a rich set covering a lot of currently proposed models of dark 
energy, on the other hand, they reveal hidden symmetries, with a particular role of the 
separatrix £l m = 2/3, and link with the issue of the content of physical constants. 
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1. Introduction 

With the establishment of observational evidence favoring dark energy domination in the 
recent history of the cosmological expansion, numerous different models were suggested 
to account for it. For most of them the value of the corresponding density parameter 
is chosen to fit the observational one. In contrast, Gurzadyan and Xue formula pQ for 
dark energy predicts the observed value for the density parameter of the dark energy. 
The formula reads jU Ej 

TT Kc 1 TT C 4 1 
8 11 ^ = g Q ~2 
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where % is the Planck constant, the Planck length is L p = (^) 2 , c is the speed of 
light, G is the gravitational constant, a is the scale factor of the Universe. Following 
the original idea by Zeldovich [Hj, it corresponds to the cosmological term 

8ttGp G x 8ttGtt c 4 1 2 f c\ 2 

Acx = ^ = -8Gr f U • (2) 

Therefore, to keep it constant, the speed of light should vary with cosmological expansion 
as c oc a. One can consider other possibilities as well, assuming the presence of some 
fundamental physical quantity and admitting variation of basic physical constants [3] 
in the spirit of Dirac approach. Since the Planck constant does not appear in (JJJ, these 
can be the gravitational constant and the speed of light. 

In the literature there is a long-standing discussion how to understand possible 
variability of physical constants, linking this point to the meaning of world constants 
and units in physics |5j. Clearly experimentally only dimensionless constants variation 
can be detected, but the meaning and the content of the underlying physical theory 
changes depending on the issue which dimensionful quantity varies. 

In fact, models with varying physical constants are among the currently discussed 
ones (e.g. |0j). However, models based on the formula (0) are very different 
from various phenomenological models of dark energy, phantoms etc., because of the 
underlying empirical basis, namely agreement between predictions for dark energy 
density parameter and the value, following from the array of recent observations. 

The classification of various cosmological models following from GX-formula is given 
in [3] . Furthermore, it is shown |7] that GX models provide indeed good fit to supernovae 
and radio galaxies data. These models contain hidden invariance, the separatrix dividing 
cosmological solutions into two classes: singular Friedmannian-type and non-singular 
with vanishing initial density [Zj,jHl- I n this paper we provide cosmological equations 
for the simplest cases discussed in j3] and generalize them, including radiation field, and 
compare them to the standard cosmological model. 

In the section 2 cosmological framework with varying physical constants is 
formulated. In the section 3 we derive cosmological equations for particular cases, 
discussed in jlj. We provide analytical solutions for energy density in all models, and 
also in some models for the scale factor. Conclusions follow in the last section. 
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2. Cosmological equations with varying constants 

In this section we are going to derive cosmological equations, assuming varying 
constants, such as the speed of light c(i), the gravitational constant G(t) and the 
cosmological term A(t) depending on the cosmic time t. Einstein equations are 
postulated to be valid in such 

We start from the FRW interval 

dr 2 



ds 2 = -c(t) 2 dt 2 + a 2 (t) 



r 2 



(d6 2 + sin 2 Odip 2 



(3) 



1 — kr 2 

where r, 9, ip are spatial coordinates, a is the scale factor. Using Einstein equations 

and the energy-momentum tensor 

T^v = diag(-e,p,p,p), 

with p being the energy density, p being the pressure, we arrive at the cosmological 
equations 



tt2 kc 2 A 8nG 
or 3 3c z 



H + H 2 -± = -^(p + 3p) + H C -, (5) 

3 3c 2 c 

where k is the sign of spatial curvature, H = dlna/dt is the Hubble parameter, a dot 
denotes time derivative. From (jlj) and J5J) (or from ^f-T^ u . u = 0) continuity equation 
can be written 



p + 3H(p + p) = -p A + {p + p A )[4---), (6) 




where have we introduced p A = Ac 2 /(8nG). The same equations were derived 
e.g. in PJ. The difference with jH] is that we consider the scaling (0) not as a 
phenomenological, but as an empirical relation. 

It is easy to check that (jlj) is the first integral of ((SJ) as it is within standard 
Friedmann cosmology. Cosmological equations for the case when c and G both depend 
on time can be found in [10]. For the case when instead A and G depend on time 
these equations are given in [TT]. Our equations reduce to equations of [TU] and [TT] 
with the corresponding assumptions. There is a little difference with [TO], because we 
postulated only invariance of the Einstein equations, while authors of JU] assumed that 
both Einstein and Friedmann equations are still valid. With these assumptions the last 
term in © disappears and our equations looks identical to those of [TU] . 

Note, that the former energy conservation does not hold any more: there are sources 
and drains of energy (e.g. from dark energy to the usual matter, from the variation of 
the speed of light and the gravitational field strength) in (JUJ). As we will see in the next 
section, (jHJ can be integrated to obtain energy density as explicit function of the scale 
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factor. Further for some models Friedmann equation (jlj) can be explicitly integrated to 
get the scale factor as a function of time. 

In j3] five different cases have been discussed, namely: 

(i) c=const, G=const, A ^const; 

(ii) coca, G=const, A =const; 

(iii) c=const, G oc a~ 2 , p\/c 2 =const, A ^const (case 3.1); 

(iv) c oc a 1 / 2 , G=const, p\ =const, A ^const (case 4.2); 

(v) c ^const, G ^const, A ^const. 

In what follows we analyze each of the first four cases, provide cosmological 
equations and corresponding solutions for them. 

3. Cosmological models based on Gurzadyan-Xue formula 

In this section we present cosmological equations for each of the cases (i-iv) and discuss 
the corresponding cosmological scenarios, providing analytical solutions, where possible. 

3.1. Models with dark energy and pressureless matter 

In this subsection we put the pressure in © and (jOJ) to zero and introduce mass density 
p = p/c 2 . Each solution can be parametrized by the usual matter density 



"m = — = OIJ 2 ^0, {<) 



Po _ 8ttGo 
p c 3Hq 

where p c is the critical density in usual Friedmannian models, subscript denotes the 
value of each quantity today. In what follows, to define density parameters we use 
normalization condition 



E n * = 1 . ( 8 ) 

i 

where % numbers components such as dark energy, radiation, usual matter. 
In order to get density as a function of the scale factor we rewrite © as 

^ + 3^ — R— * ^° (87tC/j, 4- A) p^ 1 <^A ^ 
da a A-kG cda ' Gda 8ttG da 

The crucial feature of these models is the presence of a separatrix in solutions 
for the density as a function of the scale factor (see [7] and [Hj) . This separatrix 
divides solutions in two classes: the one with Friedmannian singularity, and the one 
with vanishing density and non-zero scale factor in the beginning. 
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To see the difference from our models and the corresponding Friedmannian ones 
we compute for each case mass transfer rate, namely the ratio between R of and 
expansion term 3/i/a, i.e. 

. , O-R , . 

M = -. (10) 

If dimensionless M. vanishes, it means the usual Friedmannian solution /i oc a -3 
holds. When M. = 1 mass transfer cancels the effect of expansion and the mass density 
does not change in course of expansion. If it is large, effect of expansion can be neglected. 

Model 1 (varying cosmological constant). Neither the speed of light nor 
the gravitational constant vary with time, but in this case from (0) we find A oc a~ 2 . 
Cosmological equations are 

tt2 k'c 2 8ttG 

where k — k — ^ . It is clear that k! < 0, i.e. the effect of A variation is to make spatial 
geometry negatively curved, even if the original sign of k is " + " . The source term in 
the continuity equation can be interpreted as the energy transfer from the dark energy 
component into the matter component (see ^T] and references therein). 
Solution for the density is given by 

and after substitution of this solution to (jJJ) we obtain 

3 



where 



±A{t = ^aA 1 (aA 1 + B x ) - ^a A 1 (a A 1 + B x ) + 

|D1 ( y/^A^ + y/aoAx + BA (12) 
+Bl ' 0g [ VaA 1 + VaA 1 + Bl J' 

A\ = c 2 (tt 2 — k), B\ — — a (4G/i ao — ttc 2 ). 

3 

This solution reduces to a simple form for separatrix value with the condition B\ = 0, 
namely 

a = ao + \J~A~it. (13) 

The separatrix solution is given by the density parameter Q m w 2/3 and depends 
slightly on the spatial curvature. For high values of the density parameter solutions 
practically coincide with the corresponding Friedmannian ones. In such a case one can 
neglect the source term in the continuity equation. This also can be done for large scale 
factors since solutions approach the separatrix where they describe Friedmann Universe 
with negative curvature. 
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The present scale factor is well defined quantity even for k = models 



a ° = ^TrTT^TT' (14) 

since this is the model with effective negative curvature with = — ■ With this 
definition ((BJ) is satisfied along with with (J2J). 

Mass transfer rate (jl(Jj) for this model in terms of density parameter is given by 

M- — (— - 

and it is clear that M. gets infinite with Q m — > 0, i.e. vacuum mass transfer 
dominates expansion, and in the opposite case the model approaches Friedmannian 
matter dominated one with negative curvature. 

Model 2 (varying speed of light). To keep interpretation of A as a cosmological 

/ » \ 1/2 

constant, we have to require the varying speed of light, c = ( ^ ) a. At the same time 
the gravitational constant does not change, G=const. Cosmological equations reduce 
to the following system 



H 



2 



/i + 3Hfi 



A' _ 8nG 

A 




where A' = A (l - f§) . The source in the continuity equation now comes from the 
speed of light variability. 

There is complete analytical solution for this model 

fi(a) = fi Q — + (l - — \ (15) 



a AixG V a J 



a(t) = - A + A2+ 2 2 £° B2 cosh (tVB 2 ) ± #(A 2 + a 5 2 )sinh (tVB 2 ) . 



where 



A 2 = -oq (4vrG/i - A) , B 2 = A[l--\. 

In this model the separatrix is given by Q m = 2/3 independent on the curvature. 
If one neglects the source in the continuity equation (that is again possible for large a 
or large fi m ), which consists of two terms, proportional to H 3 and HA respectively, we 
come to the usual Friedmann-Lemaitre cosmology with the cosmological constant || A'. 

The density parameter for A term is 

n A 
• 




and the scale factor today is arbitrary quantity in this model. 

|| Note that the cosmological constant appears naturally within McCrea-Milne approach \12\ . 



Cosmological models with Gurzadyan-Xue dark energy 



7 



Mass transfer rate (|1U|) is simply 
M = —, 

so it precisely cancels expansion with Q m — > 1. 

Model 3 (varying gravitational constant). We introduce the fundamental 
constant \lgx — s^; =cons ^' ^ ne density of the vacuum. In this case the gravitational 
constant changes, G = (^j , while the speed of light does not. The cosmological 
term also acquires dependence on the scale factor through G and is given by ©. 
Cosmological equations are 

a z 3 \a/ \ fi GX ) 

fi + 3Hfi = 2Hfi (l + , 

These are very different from the usual Friedmann equations. Even if could neglect 
the source terms in the continuity equation (that obviously cannot be done), the first 
cosmological equation does not correspond to the first Friedmann one. However, phase 
portrait of solutions for this model on the plane {/z, a} looks similar to the one of model 
2 0. 

Again there is solution for density 



/i = no h 2fi GX 1 

a V a J 



Notice similarity with (|T3j) . This is due to the fact that A is a constant in the model 
2, while Hex is constant in the model 3, so the dependence of matter density on the 
scale factor is the same, although the dependence of dark energy on the scale factor is 
different. 

Solution for the scale factor looks instead exactly the same as the solution for model 
1, (HU with 

7T 2 

A 3 = Ai, B 3 = c 2 a (/i - 2hgx)- 

^Hgx 

The separatrix is the same as in model 1, given by B 3 = 0. 
Matter density parameter for this model is 

7T 2 c 2 h 

U urn 



3 HqOq hgx ' 

that is the same as The scale factor today is the same as in model 1 and is given by 
(I14|) . as this expression follows from the definition of density parameter for curvature, 
which is the same as in model 1. Another, equivalent expression follows from definitions 
of vacuum mass density and critical density and it gives 
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Mass transfer rate (|1U|) for this model is 
M 2 



30 m 

that is analogous to the model 2, but with Q m — > 1 it corresponds to a constant fraction 
of expansion rate. 

Model 4 (varying speed of light, constant vacuum energy density). The 

constant vacuum energy density is introduced as pax = |^| =cons t; requiring variation 
of the speed of light c = (^f^f^ a 1 / 2 . Here G=const and A = 2 v / 2vr 3 / 2 v ^^a- 1 = 

3/3 



—, — 3fey, where the constant (3 appears in the cosmological equations 

* = (IT) 
6 a 

with/3 = ^(G PGX ) 1 / 2 (l-f|). 

The source term in the continuity equation can be neglected for large a, but the 
first Friedmann equation contains the positive a -1 term which is again very different 
from the standard cosmology. 

There is complete solution for this model with 



and 



Ba 



where 

^4 = a o (A*oao - y ~2Q- J > B 4 = 2d — - — (tt — k). 

This solution reduces to (fT3"j) for £> 4 = or, in other words, when O m — > 1. The 
separatrix is given by O m = 2/3, like in model 2. 

Here new density parameter for dark energy is defined as 

n - p 

ti a 

and the scale factor value today follows from this definition 

C 7T 



a 



Bo ^3(1 - n m ) 

Notice that it is different from the scale factor in models 1 and 3 by the factor 
1 - 3k/n 2 . 

Mass transfer rate (|10|) for this model is given by 

30™ + V + 7T 2 ) V0 m 
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from which we find that with Vt m — > the vaccum again dominates expansion, while in 
the opposite case this is just a constant fraction of expansion rate. Therefore, it cannot 
be reduced to the usual Friedmannian model. 

Finally, in the case 5 of [I] cosmological equations are given by the system ((H)-©. 
Definitions (jZJ and (jHJ) can be used in this general case as well. 

3.2. Models with dark energy, dust and radiation 

So far we considered cosmological models based in GX formula (JTJ), containing only dust 
and dark energy. In this subsection we turn to a more general class of models, including 
into consideration radiation. 

It is known (e.g. that there is ambiguity in definition of multicomponent 

models with dark energy and/or variation of constants. There exist strong bounds 
on coupling between dark energy and radiation so we avoid this ambiguity 

assuming decoupled radiation, namely 

Pr = PrO ( — ) , Pr = ~7T) (18) 

\a / S 



where the subscript "r" refers to radiation. In standard cosmology any relativistic 
particle in thermodynamic equilibrium can be described by these equations. However, 
varying speed of light implies variation of dispersion relation e 2 = lie 2 + moc 4 , where 
e and II are particle's energy and momentum and m is its rest mass, hence influences 
equation of state for ultrarelativistic fluid. In this paper we constrain ourselves with the 
case of thermal photons only. 

Radiation density parameter in what follows is defined by 

n 87rG ° rim 
= ^TrTaPrO- (19) 

Then total energy conservation (JHJ) implies 

dp +3 P = _A_ ( 2 — - dA ^ + 

dG\ dc ( 2 °) 



da a 8nG \ cda Ada 



p r0 f a 

P+ — ( — 

c z Veto 




Gda I ^ cda 



It is clear from (J2U|) that in general matter couples to radiation always except for 
the case c = G=const. 

In this section we provide solutions for the set of cosmological models described 
by P|).(J5|). fT8|) and (f2Uj) . Generally speaking, these models do not reduce to the ones 
considered above because of special role of radiation, given by (fTH|) . 

Model 1, with radiation. This is the unique case when the radiation does not 
couple to matter density, because neither speed of light, nor the gravitational constant 
vary. 
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Cosmological equations for this model read 

ff2 W 8^G/ + P*(±y>) t (21) 



a 2 3 y c 2 \a 

Solution for the matter density is the same (fTTJ) as it is in the absence of radiation. 
Additional peculiarity of this case is the presence of separatrix, playing the same 
role as for models without radiation. The separatrix is defined by 

2 1 - Q r 

On - 

that again gives fl m ~ 2/3 for any curvature parameter. 

Implicit solution for the scale factor as a function of cosmic time is 

t _ t0 1 (J^ _ J^) _ 4 log Ml + *A + WW 



^2 V / 2i| \ii + 2a i 2 + 2ViiF(a ) 



where 



and 



= A + Aix + A 2 x 



8itG 4 



^0 — -ITYPrOOjQ 
oC 

~ 87TG / 2 7TC 2 \ 

A 2 = (vr 2 - jfe) c 2 . 
The solution for separatrix is explicit and looks rather simple 



This solution approaches (jl3j) at late times. Solution for density is represented 
at fig. (0) for selected values of Q m . Like in models without radiation solutions with 
Q m > 2/3 reach Friedmannian singularity in the past, while solutions with Q m < 2/3 
start with vanishing density. 

Both density parameter and critical density are defined by (J7J). Radiation density 
parameter is as usual (|19|). Curvature density parameters remains the same as in 
model without radiation, and the scale factor becomes a little bigger due to the change 

Vi - -> y/i-n m - n r m gap. 

Model 2, with radiation. Since the speed of light varies in this model, matter 
becomes coupled to radiation and solution gets complicated. 




Figure 1. Density of matter in model 1 in terms of critical density, depending on the 
scale factor in terms of present scale factor with VL m equal to 0.3 (firm curve), 0.667 
(thick curve) and 0.999 (dotted curve). 



Cosmological equations for this model are 

,,n A' 8nG . , 
//" J -7- ( -/Vo-ffl- 

fi + H(/i- 




Hence the solution for the density is 



/i(a) = - 
a 



A , , 47r 2 / / a 

a ) + fi a + — — p r0 1 - — 



AixG 5Aao V ^ a o 



This model does not possess separatrix in the sense discussed above, because the 
early expansion is dominated by the negative definite term (since the vacuum energy is 
positive definite). The solution for matter density in this model is represented at fig. 
©• 

Although there is no separatrix with respect to initial singularity, the value 
fl m = 2/3 again plays important role, dividing solutions in two classes, depending 
on the sign of time derivative of the density: for Q m > 2/3 the density decreases with 
time in matter dominated epoch, like in Friedmannian models, while for fl m < 2/3 the 
density increases with time. This behavior is seen from equation for the density, since 
the condition A = AnGp leads to Q m = 2/3. 

In contrast to simple GX model without radiation, the scale factor is no more 
arbitrary in this model and is given by (fT^j) with the change y/1 — Q m — > \/l — Q m — Q r . 
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Figure 2. Density of matter in model 2. Q, m equals to 0.3 (firm curve), 0.667 (thick 
curve) and 0.999 (dotted curve). 



Radiation density parameter is 



7T" p r o 



A /leal' 
and it coincides with (fTTij) . 

Model 3, with radiation. Cosmological equations now reduce to 



H 2 



27T 2 



1 / PrO 

PGX \ c 1 



(23) 



p + H(p- 2p GX ) = 2H 



PrO ( a_ 
a 



Solution for density in this model is 

p{a) = 2p GX + —{Ho ~ 2/icx) + - -o~ 1 - — , 
a 6 c z a \ V a / J 

and its dependence on the scale factor is similar to the one shown at fig. (j2J. This is 
because again with negligible contribution of radiation the time derivative of the density 
is determined by the condition p = 2p GX which leads to the same density parameter 

Q"m = 2/3. 

The scale factor is the same as in the model 2 above. Radiation density parameter 

is 

^ = 7^ PrO 

3 pcx^Hq ' 

and it coincides with (fT^j) due to definition of pox- 
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Model 4, with radiation. For this case cosmological equations read 

H 2 = — [t*+-J^ PrO - +-, 

3 \ ay 8&Pgx Vao/ / a 

a V 2G \ p GX \a J ) 
The solution for model 4 is 

and it is represented at fig. (J3J). 



7 n 




0.5 1 1.5 2 2.5 3 a o 



Figure 3. Density of matter in model 4. £l m equals to 0.3 (firm curve), 0.667 (thick 
curve) and 0.999 (dotted curve). 

There is no special value of density parameter in model 4. Solutions start with 
zero density and positive scale factor, then the density increases rapidly in radiation 
dominated stage, but entering into matter dominated stage it starts to decrease, so for 
each solution there is a maximum in the density (see fig. (J3J)). 

Here the density parameter of radiation is defined by 

_ 16tt 3 / 2 [G~ 
ZHfiao V Pgx 
and it is again equivalent to (fTUj) . 



(24) 
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4. Discussion and conclusions 

In this paper we explored various cosmological models with Gurzadyan-Xue dark energy 
formula (JTJ, assuming also variation of physical constants such as the speed of light and 
the gravitational constant. It is shown that unlike standard Friedmann models, various 
interesting possibilities appear. In particular in the model 1 the continuity equation is 
reduced to the energy transfer from dark energy component into the usual matter. For 
some cases (models 1 and 4), the source term in the continuity equation can be neglected 
for sufficiently large scale factors (or sufficiently large present density of matter), and 
the discussion of the dynamics in this asymptotic cases can be given in terms of the 
standard cosmology. 

Solutions for the density in terms of the scale factor for all models with pressureless 
matter and dark energy contain separatrix which divides solutions in two classes. The 
first class is analogous to the standard Friedmann solutions with singularity in the 
beginning of expansion. The second class, in contrast, does not contain singularity; 
solutions start with vanishing density and non-zero scale factor. 

We also considered generalization of GX models, including radiation, assuming that 
dynamics of the latter is not influenced by the variation of constants and by dark energy. 
Interestingly, for the model 1 we obtained the behavior qualitatively similar to models 
without radiation. For the rest of models separatrix in the sense discussed above is 
absent and all solutions belong to the second class, namely non-Friedmannian solutions. 
The density in these solutions may be larger in the past, but if so, it has a maximum 
and all solutions start with vanishing density in radiation-dominated epoch. However 
for models 2 and 3 with radiation the same value of density parameter fl m = 2/3 that 
parametrizes separatrix in simple GX models, corresponds to constant mass density of 
pressureless matter in matter dominated epoch. For larger density parameter the mass 
density decreases with expansion and for smaller density parameter it increases. 

In fact, not all GX models with non-Friedmannian behavior (Q m < 2/3) pass 
age constraint and have deceleration parameter which fits the value followed from 
observations. At the same time model IV with Q m > 2/3 pass both constraints and 
represent a good candidate for dark energy model j7j which in the limit of vanishing 
energy transfer from dark energy to the usual matter has a~ l term in Friedmann 
equation. In the same way models I and III fit these constraints for density parameter 
slightly smaller than the separatrix value. This circumstance once more points out on 
a special role of the separatrix in GX models. 

Notice another interesting feature of GX models with radiation. It is well known 
that Friedmannian models with zero curvature have arbitrary scale factor. The same 
holds also for simple GX models (model 2). However, this degeneracy breaks down in 
GX models with radiation, where the scale factor can always be computed in terms of 
present density parameters, and present values of (varying) physical constants. 

We are thankful to the referees for helpful comments. 
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